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Abstract. The total single-photon ionisation cross section was calculated for helium 
atoms in their ground state. Using a full configuration-interaction approach the 
photoionisation cross section was extracted from the complex-scaled resolvent. In 
the energy range from ionisation threshold to 59 eV our results agree with an earlier 
_B-spline based calculation in which the continuum is box discretised within a relative 
error of 0.01% in the non-resonant part of the spectrum. Above the He ++ threshold 
our results agree on the other hand very well to a recent Floquet calculation. Thus our 
calculation confirms the previously reported deviations from the experimental reference 
data outside the claimed error estimate. In order to extend the calculated spectrum 
to very high energies, an analytical hydrogenic-type model tail is introduced that 
should become asymptotically exact for infinite photon energies. Its universality is 
investigated considering also H~ , Li + , and HeH + . With the aid of the tail corrections 
to the dipole approximation are estimated. 
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1. Introduction 

Since the beginning of quantum mechanics the photoionisation cross section (PCS) 
of the helium atom was investigated in a number of experiments and numerical 
calculations. Helium is one of the simplest quantum mechanical systems that is relatively 
easily experimentally accessible, but also amenable to very accurate calculations. 
At the same time, it is theoretically challenging, since even within non-relativistic 
quantum mechanics the helium atom cannot be solved analytically. Due to this special 
characteristics the helium PCS is very attractive for a direct comparison of theory and 
experiment. 

In 1994 high-precision measurements of the PCS were performed by Samson et 
al [U [2] using a double ion chamber and a high-voltage spark discharge. Since then 
the therein reported values for the PCS of helium with an estimated accuracy of 
1 — 1.5% in the low-energy range and about 2% for the high-energy part beyond 
the double ionisation threshold were used in many applications in, e.g., astrophysics, 
plasma physics, and chemistry. Recently, the PCS of helium became also relevant 
for the characterisation of novel light sources like high-harmonic radiation or free- 
electron lasers (FEL). The supposedly very accurately known PCS of helium provides 
a natural way for tests and calibrations (especially of the intensity) of these new- 
generation light sources [3j H]. For example, in the SASE (Self- Amplified Stimulated 
Emission) experiment at FLASH (Hamburg) the two-photon double photoionisation 
of helium is used to determine the duration of ultrashort femtosecond pulses [5]. 
Thereby, the nonlinear autocorrelation of direct He — > He 2+ + 2e~ and sequential 
He — > He + + e~ — )■ He 2+ + 2e~ double ionisation process is measured. The first step of 
the sequential process corresponds, of course, to single-photon ionisation of He. 

However, despite the development of new theoretical approaches and the access 
to increasing computational power it has so far not been possible to reproduce the 
experimental reference data in [TJ [2] to the therein claimed accuracy. Supposedly 
very accurate theoretical values for the PCS of helium from the ionisation threshold 
to photon energies of 71 eV were reported by Venuti et al in [B] . This was a follow- 
up work to the one presented in [7] which extended to 2keV but was less accurate. 
The calculations were performed within the configuration-interaction approach in which 
the orbitals were expressed in B splines (for the radial part) multiplied by spherical 
harmonics (for the angular part). The finite range of the adopted B splines leads to a 
box-type discretisation for the continuum wave functions. From convergence studies and 
the agreement between the results obtained using the length, velocity, or acceleration 
form of the dipole operator the authors of [6] estimated the error to be smaller than 
0.001 Mb which corresponds to a relative error of 0.014 — 0.063 %. In comparison to the 
experimental data in [U |2] a deviation of up to about 2.6 % was, however, found in the 
considered low-energy range which is almost twice the estimated experimental error. 

More recently, Ivanov and Kheifets implemented a Floquet approach and calculated 
the PCS of helium starting at a photon energy of 80 eV [8]. Again, the authors claim to 
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reach an accuracy of the order of the fraction of a percent, but find noticeable deviations 
from the experimental reference data that easily reach up to 6%. For a single value of 
the photon energy, 40 eV, a good agreement was on the other hand found with the 
theoretical results in [6]. As a consequence, the non-relativistic dipole and infinite- 
mass approximations used in both calculations appear to be inadequate or either the 
calculations or the experimental data are less accurate than claimed by the respective 
authors. In order to shed more light on this question, we performed calculations with a 
different theoretical approach in which the orbitals used in the subsequent configuration- 
interaction method are constructed from Slater-type orbitals. Furthermore, the PCS is 
extracted from the complex-scaled resolvent and thus no box discretisation as, e.g., in 
[6] is used. We extend our ab initio results to very high photon energies by introducing 
an analytical, hydrogen-like model tail. This is of interest, for example, in view of the 
present efforts of extending the FEL sources to the x-ray regime, as with the LCLS (linac 
coherent light source) at Stanford, the XFEL (x-ray free electron laser) in Hamburg, 
or the SPring8-XFEL project in Japan. The universality of the here introduced tail is 
investigated by considering also the other two-electron systems H~, Li + , and even the 
molecular ion HeH + . The tail is finally used to obtain an analytical estimate of the 
first-order correction to the dipole approximation. 

2. Method and Computational Details 

2.1. Ab initio calculation 

In view of the large mass difference of the He nucleus and the electrons we adopt the 
infinite-mass approximation for the He nucleus in our calculations. This should be 
justified, since for the PCS of the He + ion, for which the analytical result is known, 
the size of the modification due to the finite mass of the nucleus is only of the order 
of 0.1%. It is expected that this effect has a similar size in the case of the neutral He 
atom. Furthermore, there should not be a strong energy- dependent contribution due to 
the finite mass. In fact, for the ion it is energy independent. Relativistic effects should 
also be negligible, because their scaling parameter, the nuclear charge, is evidently a 
small value in the case of helium atoms. 

To obtain the wavefunctions and the corresponding energy eigenvalues of He we used 
a direct expansion in Slater-type orbitals (STOs). In the calculation of the X S helium 
ground state the same basis-set parameters were used as in our recent calculation of the 
final-state spectrum of helium atoms after /3 decay of tritium anions [9]. The previous 
results for a large number of energy eigenvalues exhibited a very good agreement with 
literature values. For the 1 P states we modified the X S basis set as to adapt it to the 
different symmetry. As a result, we used 360 STOs leading to 3331 configuration state 
functions (CSFs) for all *P states and 555 STOs resulting in 3481 CSFs for the X S 
ground state. More details about the structure of the basis set and the construction 
of the symmetry-adapted CSFs used in the subsequent configuration-interaction (CI) 
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calculation can be found in [9]. 

The photoionisation cross section a is related to the optical oscillator strength 
density df/dE by the relation [TO] 

ne 2 h df(E) , nnncnniin .df(E) 



a 



109.7609 Mb -^eV (1) 

2e m e c dE dE 

with the electron mass m e and charge e, the reduced Planck constant h, the electric 
constant e , and the speed of light in vacuum c. Using atomic units (m e = 1, e = 1, 
h — 1) the df/dE can be evaluated from the dipole transition probabilities P{E) 

( df(E)\ f EP(E) in length form 

\ dE ) a u ~ j p(E) jE in velocity form. 

Finally the P{E) can be extracted from the complex-scaled resolvent according to 

mm 

(*un\d (mm)(n(nd( mmy 

El{6)-Ef{e)-E 

The complex-scaling angle is denoted by 9 and (ty(9*)\ is the biorthonormal eigenstate 
to \^/(9)). It is obtained from the latter by a transposition and complex conjugation 
of the angular part, while the radial part is only transposed but not conjugated. A 
variation of the angle 9 allows the determination of an optimal angle 6^ by requiring 

9P(E) . 

- mm. (4) 




89 

In other words, the best approximation of P{E) is obtained with that value of 9 for which 
P{E) shows the smallest dependence on 9. For more details about complex-scaling see 
[HI HU H3] and references therein. An application of this method for the photoionisation 
of the molecular ion HeH + was described in |12j . 

In Equation the are the final helium l P wavefunctions and E^ the 

associated energy eigenvalues. Analogous definitions apply to the X S initial state. The 
operator d describes the coupling of an atomic electron to a (classical) electromagnetic 
field. Thereby the latter is represented by a plane wave in the spatial domain, e^ r , 
with the wave vector k and the spatial vector r. Without loss of generality we can 
assume the wave is propagating along the z axis and the operator becomes in length 
and velocity forms 

dj = -e lk ' z z (5) 

and 

d„ = -e lk H z , (6) 

respectively. For sufficiently low energies the dipole approximation e lkz ~ 1 can be 
applied to a high level of reliability. The complex-scaled versions of the operators are 
simply given by d/(#) = —e +w z and d v (9) = —e~ w v z . For an iV-electron system the 
electronic part of the operator is simply the sum of the iV one-electron operators. 
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2.2. Analytical model for high photon energies 

For very high photon energies we introduce an analytic tail. The concept behind its 
construction is that in the case of single-photon ionisation the ejected electron takes 
away a large fraction of the energy of the absorbed photon. Thus it escapes very fast 
from the nucleus and the remaining (spectator) electron. Therefore, the fast electron 
experiences the remaining system (with a maximum screening of the nucleus by the 
spectator electron) to a good approximation as a point charge Zf which is the sum of 
the charges of all remaining particles. (In the case of a He atom, the escaping electron 
experiences the remaining He + ion as a point charge with Zf = +1.) As a result, we 
may approximate the X P two-electron wave function as (the spin part is omitted for 
better readability) 

\& sS P(g)} = 2"3 [\ls z *£ P z f) + \£p z fls z ')] . (7) 

Here, \Sp Zf ) is the Coulomb continuum p-wave function for energy £ and charge Zf, 
while \ls Zt ) is a hydrogen-like ground-state wave function with effective charge Zi. The 
initial 1 S state is approximated as a product of two hydrogenic s orbitals with the same 
effective charge Zj, 

|^ ls2 > = \ls z <ls z >). (8) 

In this model it is assumed that during the (fast) photoionisation process (escape 
of the emitted electron) the spectator electron has no time to relax. As a consequence, 
the effective charge Zi of the spectator electron remains the same in the initial flE]) and 
the final ([7]) wave functions, but changes from Zi to Zf for the fast escaping electron. 

Using this model we derive an analytical formula for the photoionisation cross 
section for two electron systems at high photon energies which reads in atomic units 



512 exp (_ ^/arctan(«/Z0 \ E ft z (V _ 2Z f (ft + 2\ 



3(l-exp(-^))(^ + ^) e 



with k = J2(E — I), the photon energy E, the ionisation potential /, and the number 
of electrons N. This expression is obtained when using the length-form representation of 
the dipole operator. For the velocity form of the dipole operator the model cross-section 
is 




*° (g > = U(V-2v ' (10) 

For the one-electron case (N = 1, Zi = Z v = Zf) the two representations become equal. 
This is, however, not the case for our two-electron model, since it adopts approximate 
wave functions. As a consequence, the corresponding commutator relation leading to 
equivalence between length and velocity forms of the dipole operator is not fulfilled. 

Within the model, two of the three parameters entering Equation (Q and Equation 
fTlOl) . Zf and /, are fixed by the sum of the charges of the remaining particles and 
the either experimentally or theoretically known ionisation potential, respectively. The 
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choice of Z^ (Z v ) is on the other hand less evident. One possibility is to relate it 
to the ionisation potential via Z{ = a/2J which leads to a parameter-free model tail 
[HJ [15]. A second choice is based on the requirement to yield the best possible 
ground-state description when adopting the trial wavefunction of © in a variational 
calculation. For He this results in the mean-field value Zi = 1.6875 [16] compared 
to about 1.34 for Zi = \p2J. A third alternative for the choice of Zi is to fit this 
parameter in order to provide the best agreement of the resulting tail with some ab 
initio or experimental PCS, within some energy interval. It should be noted that for 
N = 1 and Zi = Zf = \/2I Equation (jUJ) is also evidently applicable for H-like systems 
and in this case one obtains the result given in [T5] and proposed as an approximate tail 
also for the generalised oscillator strength density in [14J. Clearly, as is discussed below 
for the example of HeH + , also the present tail can be adopted to molecular systems. 
In the spirit of the Born-Oppenheimer approximation it may be useful to use in such a 
case an internuclear-separation dependent electron binding energy I(R) instead of the 
experimental ionisation potential [T7J [TBI [T9] . 

Through its simplicity this model offers the possibility to study analytically effects 
beyond the dipole approximation. Since the model becomes better for higher energy it 
is valid in the regime where the wavelength of the electric field may approach or even 
extends below the magnitude of 1 nm and is thus comparable to inner-atomic distances. 
To estimate these effects in helium we derived expressions for the first-order corrections 
to the dipole approximation. In this case the one-electron interaction operator in length 
form is given by 

d; = — zexp(ifcz) £3 — z — iEai 2 (11) 

with the photon energy E = chk = (k/a) (in a. u.) and the fine-structure constant 
a. Since the initial state given in Equation (jSJ) has X S symmetry, the correction term 
lEaz 2 couples it only to final states with either X S or X D symmetry. In the spirit of the 
model one has then to substitute \£p Zf ) in Equation (J7J) with \£s Zf ) and \Sd Zf ) for the 
continuum states with X S and X D symmetry, respectively. 

Since two distinct final states are reached at a given energy, the (incoherent) 
superposition of the cross-sections into the *S and *D channels is given by 

a^E) = <Jq(E) (l + (Ea) 2 (A 2 S + A$) (12) 

with 



3Zi (2Zj -K 2 -Zi (3Z f + Zi)) + Z f (2K 2 - Zj) 
^ ~ V3(Z f -2Z i ) yJZ 2 f + tf(Z? + tf) (13) 

and 

4 (Z f - 3Zi) yJZj + Ak 2 

r 15{Z f -2Z l ) (Z 2 + k 2 

In the case of the X S final states one obtains formally a term (£s Zf \ls Zi ) (ls z% \z 2 \ls Zi ) 
which can be interpreted as a shake up of the electrons. It is an artifact from the unequal 



Ad = -^-^ • (14) 
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treatment of the two electrons in our model. Since in this case only one bound electron 
interacts with the electromagnetic field, the other electron (which will be ejected) plays 
the role of the spectator and consequently experiences no relaxation. Thus to stay 
consistent within our model we have to exchange the effective charges in this term to 
(Ss Zi \ls Zi )(ls Zf \z 2 \ls Zl ) which due to orthogonality of the \Es Zi ) and \ls Zi ) orbitals 
vanishes. 

We successfully tested the derived expressions for the cross-sections in the one- 
electron case Z = Zi = Zf with the aid of the sum rule [20] 

So = (E n — E ) 

n 

In Equation ( fl5l) the sum runs over all bound states \n) and |0) denotes the initial state, 
e.g., the ground state. In the case k — 1 (dipole approximation) one has a L = <7q and 
the sum-rule result should be equal to the number of electrons. For k = 2 (corrections) 
one finds S = 4/Z 2 by using a L = a^(A 2 s + A 2 d ). 

3. Results 

3.1. Low-Energy Photoionisation Cross-Section 

In Table [1] the calculated ah initio PCS values are listed for the length and velocity 
forms of the dipole operator. A remarkable agreement of the order of 10~ 3 Mb can be 
noticed between these two formulations that for exact wave functions yield identical 
results. The results are given for the optimal values 6 opt obtained according to Equation 
(jlj). However, overall, only a very small variation of the PCS with 9 is found which 
indicates that the used basis set is rather complete for the considered energy range. 

Table [1] also compares the present results with the experimental data of Samson 
et al [lj and the supposedly most accurate previous calculation of Venuti et al [6J. If 
shown graphically, as in Figure HJ there are almost no visible differences between the 
various theoretical and experimental results. 

The very good agreement of the present results with the ones in [6] is confirmed 
in Figure |2] that shows the relative deviation to the present length-form data for 
photon energies below 59 eV. Despite the completely different theoretical approaches 
the agreement is remarkably good, since the deviation is less than 0.16% and especially 
in the energy range from 38 eV to 58 eV the relative deviation is of the order of ±0.03% 
(except for the values in [6] at 42 eV (length form) and at 52 eV (velocity form), which 
we suggest to be typos). Figure |2] shows also the deviation between our results in 
length and velocity form. Again good agreement with deviations less than ±0.04% is 
found. In an earlier work, Chang and Fang adopted in (21] a theoretical approach that is 
practically identical to the one used in [6] , but a smaller basis set was used. Our results 
agree also well with the ones in [21] (relative error of about 0.1%), but are consistently 
in better agreement with the ones in [6]. Clearly, on this level of accuracy convergence 
of the adopted basis set is finally decisive, while the two completely different approaches 



(0|z fc |n> +/ a L (E)dE = k 2 (0\z 2k - 2 \0) 



(15) 



Accurate photoionisation cross section for He 



8 



Table 1. Total photoionisation cross sections a (in Mb) for He as a function of the 
photon energy E (in eV). The present results are compared with experimental values 
of Samson et al [T] and with the B-splines calculations of Venuti et al [6] . 

This work Venuti et al [6] Samson et al [1] 
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Figure 1. (Colour online) Total photoionisation cross section of helium ( 1 S ground 
state) from the ionisation threshold to 59 eV (red — — : present work, x : theoretical 
values of Vcnuti et al |6], □: experimental values of Samson et al [TJ [2])- The inset 
shows the theoretical PCS of this work and experimental values of Samson et al on a 
logarithmic scale for higher photon energies. 



(B splines with box discretisation vs. Slater-type orbitals with complex scaling) appear 
to yield identical results within the achieved level of convergence. This indicates the 
correct and numerically stable implementation of both approaches. In fact, in the case 
of the single data point (at 40 eV) given within the energy range shown in Figure |2j the 
result of the Floquet calculation by Ivanov and Kheifets [8] agrees also within 0.03% 
with our results and the one in j6]. Therefore, three different theoretical approaches 
agree within an extremely small relative error. 

On the other hand, the experimental helium PCS of Samson et al [U [2] is still the 
primary reference for benchmarks and comparison, despite the fact that recent numerical 
calculations [6j |8] yielded cross sections that differ from the experimental values outside 
the estimated experimental error bars. In view of the already discussed very good 
agreement with the theoretical data in [6] (Figure [2]), we confirm the deviations to the 
experimental data in [TJ [2] outside the error bars. This can clearly be seen in Figure 
|3j For photon energies between 51 eV and 55 eV the relative deviation of the results of 
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Figure 2. (Colour online) Relative deviation of our results in velocity form (red □) 
and of the results of Venuti et al [6] (<)) in length (orange), velocity (black), and 
acceleration forms (blue) from the helium photoionisation cross section calculated in 
this work within the length formulation. 



Samson et al from our data and also from the ones of Venuti et al is evidently above the 
error of ±1.5% estimated by Samson et al [HE]- The largest deviation of approximately 
2.6% in comparison to our values occurs at 53 eV and 54 eV. These deviations are 
too large to be explained by a failure of the approximations adopted in the present 
calculation. Neither effects due to the finite size of the nucleus nor relativistic effects 
should be of a magnitude that is sufficient for explaining such a discrepancy that in 
addition would have to be strongly photon-energy dependent. As is discussed below on 
the basis of the derived analytical high-energy tail, also the consideration of non-dipole 
terms and thus a break-down of the dipole approximation yields corrections that are 
orders of magnitudes smaller than the ones required to find agreement between theory 
and experiment. 

In view of the in the Introduction discussed importance of the He PCS there exist 
also data sets that represent a compilation of experimental and theoretical data and 
try to cover large photon-energy ranges. Such a data set was reported by Yan et al 
[22] and was claimed by its authors to be reliable for all energies. The comparison in 
Figure |3] shows that agreement with the here considered theoretical results is reasonable, 
but not really good. In fact, at both ends of the shown energy range the agreement 
of the theoretical data with the compiled ones is less good than the one found for 
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Figure 3. (Colour online) Relative deviation from the helium photoionisation cross 
section calculated in this work (length form): theoretical data of Venuti et al [6] (orange 
<C>, length form), Chang and Fang [21] (yellow □, length form), as well as Ivanov and 
Khcifets (single point at 40 eV, blue A) [8], the experimental values of Samson et al 
[HE] (violet O ), and the compiled data of Yan et al [22] (red *). The shaded area 
illustrates the error range estimated by Samson et al for their experiment (±1% until 
48 eV and ±1.5% starting from 49 eV). 

the experimental data in [TJ |2]. Since the compiled data lie below the theoretical 
ones for low energies and above for larger energies, the agreement is only good for 
intermediate energies close to the crossing point at about 42 eV. Especially close to the 
ionisation threshold the experimental data in [U |2] appear to be clearly superior to the 
compiled ones in [22]. In fact, within the first 10 eV above the ionisation threshold the 
experimental data are in remarkable agreement to theory with a deviation of less than 
about 0.15%. 

3.2. High-Energy Photoionisation Cross-Section 

Motivated by the recent work of Ivanov and Kheifets [8] who claim an even larger 
deviation from the experimental data of Samson et al for larger photon energies than 
for the lower ones considered in [6], we also studied the single-photon ionisation process 
in the non-resonant energy regime above the He + ionisation threshold (E 79 eV). 

Although our approach, the expansion of the two-electron wave function in Slater- 
type orbitals, was originally developed for bound-state transitions, the extension by the 
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Table 2. Total photoionisation cross sections a for He at high photon energies E (in 
eV). Until E — 250 eV the a values are given in Mb and our results stem from the full 
ab-initio calculation. Starting from 400 eV the results for the model tail (|9]) are given 
and all a values are in Barn. Our results are compared to the experimental values of 
Samson et al pQ , the B-spline calculation of Decleva et al if] , and the Floquet results 
of Ivanov and Kheifets [8] . 



E 


This work 


Decleva et al [7] 


Sftirmnn pf nl 1 
Ocxp 


Tvanov pi nl IR1 
0F1 










80 


73791 


73777 


75Q 


74 


0.693 


0.7369 


85 


63041 


63099 






0.595 


0.6308 


91 


59708 


52762 






0.502 


0.5272 


95 


0.47011 


47063 






0.45 


0.4701 


100 


40960 


41 004 


41 7 


403 


0.393 




111 


30808 


30854 






0.3 


0.3082 


120 


24809 


24863 


0.251 


0.243 


0.244 




140 


0.16041 


0.16108 


0.167 


0.158 


0.160 




160 


0.10924 


0.10985 


0.112 


0.108 


0.108 




180 


0.07749 


0.07799 


0.802 


0.077 


0.076 




205 


0.05270 


0.05321 






0.051 


0.0529 


250 


0.02881 


0.02939 


0.0306 


0.0293 


0.0277 




400 


6561 


6812 


7270 


7009 


6370 




600 


1785 


1871 


2001 


1940 


1770 




1000 


335 


354 


395 


384 


339 




2000 


33.1 


35.2 


39.5 


40 


34.8 




3000 


8.42 


8.96 






8.77 




4000 


3.17 


3.38 






3.20 




5000 


1.48 


1.58 






1.47 




6000 


0.79 


0.85 






0.77 




7000 


0.47 


0.50 






0.45 




8000 


0.30 


0.32 






0.28 





complex-scaling method provides also an excellent description of high-energy continuum 
states of the He atom [9]. Figure H] shows the deviation of various theoretical, 
experimental, and compiled data from our results for energies between 80 eV and 205 eV. 
The agreement with an older B-spline calculation by Decleva et al [7j is by far not as 
good as the one found with the later work [6] of the same authors that concentrated on 
the low-energy regime. However, the deviation of our values from the Floquet results 
of Ivanov and Kheifets [8] is for most of the data points less or equal to about 0.03%. 
Only at the highest energy, 205 eV, a deviation of 0.37% is found. Therefore, we again 
confirm the discrepancy between previous theoretical calculations and the experimental 
reference data of Samson et al, this time reaching to about 6% at 80 eV, as can also be 
seen from the comparison of the different results in Table 2. 

On the other hand, we find in the energy range between about 110 eV and 160 eV a 
deviation from experiment that remains basically below 2 %. Furthermore, the compiled 
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80 100 150 200 250 

Photon energy [eV] 

Figure 4. (Colour online) Relative deviation from the helium photoionisation cross 
section calculated in this work: Floquet calculation of Ivanov and Kheifets [5] (blue 
A), B-spline calculation of Decleva et al [7] (orange <0) experimental data measured by 
Samson et al p] (violet O ) or Bizau and Wuilleumier [23] (green <£>), and the compiled 
values by Yan et al [23] (red *). The shaded area illustrates the error range estimated 
by Samson et al for their experiment (±2%). 

data of Yan et al [22] agree in this complete energy interval better with the theoretical 
results than the experimental results of Samson et al, in contrast to the findings at 
lower photon energies. In the energy interval between about 80 eV and 100 eV the 
experimental values reported by Bizau and Wuilleumier [23] deviate from our calculation 
in a qualitatively very similar fashion as the experimental data of Samson et al. However, 
quantitatively, the deviation is smaller and remains in between 3 and 5 %. In between 
100 eV and 150 eV the data in [23] are on the other hand substantially smaller than 
our theoretical results, leading to a deviation of up to 7%. For the higher energies 
shown in Figure H] the deviation of both experimental data sets from our results is again 
qualitatively similar, but the ones of Bizau and Wuilleumier lie below our results and 
approach the latter with increasing energy, while the ones of Samson et al lie above and 
thus agreement with our data becomes worse for increasing photon energy. 

For very high photon energies above 300 eV we propose to use the analytical 
model tail introduced in Equation (|9]). In order to determine the tail parameter Zj, 
we performed a least-squares fit of the tail to our ab initio calculation in the energy 
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range from 136 eV to 272 eV and obtained Z± = 1.5293. This value lies in between 
the prediction based on either Zj = \/2I = 1.34 or the mean-field result, Zj = 1.6875. 
Although the fit interval spans only a small energy region of the calculated PCS, the 
extrapolated PCS agrees in a much larger energy range well with the ab-initio results and 
in the complete shown energy range from 300 eV to 8 keV very well with the experimental 
data of Samson et al (Figure [5]). In fact, we found almost the same effective initial 
charge, Z^ = 1.53 ± 0.01, when fitting the tail to the experimental data of Samson et 
al. The indicated uncertainty arises from different starting points of the energy span 
at the fitting procedure, while the end point was always chosen at an energy of 8 keV. 
Figure [6] shows two alternative tails besides the already discussed one. If Zi is fixed 
to its mean-field value (1.6875) and Zf is used as a fit parameter, the resulting tail 
lies above the experimental data of Samson et al. While the slope differs for photon 
energies at about 300 eV, an almost constant off-set is found between the tail and the 
experimental data for large energies, if plotted on a doubly logarithmic scale. A further 
off-set is observed, if all tail parameters are fixed on the basis of simple arguments, i. e., 
the initial charge Zj is set to the mean-field value 1.6875 and the final charge Zf to the 
value 1.0 that should be appropriate at very large separation of the ejected electron. 
Not shown is the tail that is obtained under the assumption of no relaxation (Z, = Zf), 
since this tail does not agree to the experimental data at all. 

Alternatively to the described tail, an asymptotic representation of the cross section 
is often used as an estimate in the high-energy regime. The asymptotic limit of our 
length-form model tail is (in atomic units) 

a^(E) » N b^M^l E -V2 (16) 

and for the photoionisation of helium we obtain a^ e (E) rs 487.956 E (keV) -7 / 2 . The pre- 
factor is consistent with the values discussed by Samson et al for different experimental 
data. It should be noted, however, that the analytic tail in ([9]) possesses a much larger 
validity regime than the asymptotic expression in (fTB"]) . since it shows good agreement 
with experiment and ab initio theory already for much lower photon energies. 

As mentioned above the length and velocity representations of the dipole operator 
lead to different expressions for the model tail. However, one may use the different 
representations to perform a consistency check of our model and to learn about the 
reliability of the fit procedure. For energies near infinity where the model is surely 
applicable the representations should lead to equivalent values. Thus we demanded 
the identity of the length and velocity model tails in the first order of the asymptotic 
expansion. The asymptotic limit of the model tail in velocity representation is 

a^E) - N 8 -^E~^ . (17) 

Combining ( I16p and (1T7|) allows to define a fixed relation between Zj and Z v and thus 
to define a new parameter Z* that fulfils this relation, 

Z: = (4Zf-4Z / Z i 4 + Z / 2 Z ) 3 ) 1/5 . (18) 
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Figure 5. (Colour online) Analytical model tails for different representations of the 
dipole operator: length-form (red line) and velocity-form (blue line) for their optimal 
parameters Zi and Z v . Also the present ab initio results and the experimental data of 
Samson et al are plotted. The inset shows the curves on a linear scale for low energies. 



For Zi = 1.5293 and Zf — 1 this expression leads to Z* = 1.7224, while by applying the 
same fit procedure for Oq as was used before for ctq we find Z v = 1.7455. Clearly, the 
differently obtained values for Z v and Z* are in quite reasonable agreement. Thus we 
conclude that the model tail, although it is derived with approximate wave functions, 
yields for very high energies the required independence of the chosen representation 
(length or velocity form) of the dipole operator. Furthermore, because the results of 
Equation ffTB"]) and the fit are in very good agreement, it is evident that the fit procedure 
is most suitable to derive the parameter values for the model tail. It may finally be 
observed that Z v (and Z*) are close to the mean-field value 1.6875. This may indicate 
that, if one would like to avoid the fit or if no data are available for performing a fit, the 
best choice for a parameter-free tail is the tail in velocity form adopting the mean-field 
prediction for Z v . 

In Figure [5] the length and velocity model tails are plotted for their optimal values. 
It can be observed that the photon energy where the model tail approaches the ab 
initio calculation is lower for the velocity form than for the length form. In fact, it is 
quite surprising how well the velocity-form tail agrees to the full ab-initio results and 



Accurate photoionisation cross section for He 16 



10" f 




300 500 1000 2000 5000 8000 

Photon energy [eV] 

Figure 6. (Colour online) The analytic tail (see Equation ([9]), solid lines) that is an 
approximation to the photoionisation cross section of He at high energies is shown 
for different parameters Zi and Zf (as specified in the graph) and is compared to the 
experimental values of Samson et al [T] (violet O ) and the calculation of Decleva et 
al [7] (orange <)). The dotted and dashed lines show the first-order corrections to 
the dipole approximation according to Equation (TT21) due to final states with S or D 
symmetry, respectively. The inset shows the high-energy part of the spectrum on a 
linear scale. 



the experimental data even down to the ionisation threshold. This superiority of the 
velocity form compared to the length form found at low energies is, however, on the 
first glance a little bit surprising, since the transition dipole matrix elements in length 
form are usually supposed to be preferable at low photon energies [HI El] • This is often 
explained by the lower sensitivity of the length-form matrix elements to errors in the 
long-range part of the wavefunctions and the fact that low-energy transitions are more 
sensitive to this long-range part. However, it may be remembered that the obtained 
value of Z v is quite close to the mean-field prediction. Since the low-energy part of 
the photoionisation spectrum is more sensitive to the details of the atomic potential, a 
model like the velocity-form tail that is closer to the mean-field prediction may thus be 
favourable compared to the length-form tail with a rather different value found for the 
corresponding parameter 
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Table 3. Parameters for the model tail for various one- and two-electron systems. 



Parameters 


H" 


H 


He 


He + 


Li+ 


HeH + 


Zi 


0.5586 


1 


1.5293 


2 


2.5191 


1.9414 


Z v 


0.7662 


1 


1.7455 


2 


2.7390 


1.9026 


Z* 


0.7371 


1 


1.7224 


2 


2.7152 


1.9177 


Zf 





1 


1 


2 


2 


2 


N 


2 


1 


2 


1 


2 


2 


J[eV] 


0.75436 


13.606 


24.5912 


54.4234 


75.64 


45 



In order to demonstrate the generality of the analytical tail, we investigated 
its applicability to other two-electron systems. Besides the neutral helium atom we 
considered the hydrogen anion, the lithium cation, and even the HeH + molecular cation. 
The used tail parameters are listed in Table |3j Figure [7] shows a comparison of the 
tails for these systems with corresponding literature data and, for the atoms, also with 
ab-initio results calculated in this work. The found agreement is in all cases good 
or even very good. This indicates the universality of the tail concept. It is, in fact, 
interesting that the tail works well even for a system like H~ in which the escaping 
electron experiences within the tail model no influence from the remaining hydrogen 
atom, since polarisation effects are ignored. Furthermore, the initial state of H~ is only 
bound due to correlation and thus a mean-field model is in principle not applicable. 
Nevertheless, also in this case the tail seems to work well. In the case of HeH + a high- 
energy tail had been proposed before [12], but a paramter-free tail (with Zi = Zf = \/2I) 
had been used. The result is also shown in Figure It also compares reasonably well 
with the full ab-initio results. The reason is that in this case Zi is very close to Zf and, 
in fact, also to vT- Note, for a molecular system there is the additional complication 
due to nuclear motion. However, in the spirit of the high-energy tail it should have a 
negligible influence on the high-energy part of the photoionisation spectrum. Therefore, 
the tail is formally obtained for a single internuclear separation, usually the equilibrium 
distance. Here, we chose in agreement to [12] the ionisation energy / = 45 eV. 

According to (fl2|) the model tail provides the possibility to give an analytical 
estimate of the first-order correction to the dipole approximation. The two contributions 
from final states with S and D symmetry are also shown in Figure El While the D 
contribution is clearly dominant, the total correction to the dipole approximation is 
still negligibly small at photon energies of several keV. Note, however, that the relative 
contribution reaches already 11% at 8keV due to the very small cross section of less 
than 1 barn. This allows to conclude that very accurate studies of the photoionisation 
cross section of He performed in the keV photon energy range must take corrections to 
the dipole approximation into account. However, since the correction is proportional to 
the squared photon energy and the cross section itself increases by almost five orders of 
magnitude when going from 8 keV to 300 eV photons, corrections due to a break-down of 
the dipole approximation are negligible within the here considered level of accuracy for 
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Figure 7. (Colour online) Analytical model tails (length form, red solid lines) for 
different two-electron systems are compared to literature data, a) H~: Venuti and 
Decleva [25] (blue O ), b) He: Samson et al !l] (violet O ), c ) Li + : Verner et al [26] 
(blue line), and d) HeH + (parallel contribution): Saenz [12] (blue line). For the atoms 
(a to c) also ab-initio cross sections obtained within this work are shown (black solid 
lines), while for HeH + the alternative tail proposed in [T2] is additionally given (purple 
solid line). 

photon energies below 200 eV. Therefore, the deviations between the theoretical results 
of the present work as well as the ones in [25] or [8] from the experimental reference data 
of Samson et al [TJ [2] or Bizau and Wuilleumier [23] cannot be explained by a failure of 
the dipole approximation. 

4. Summary 

An ah initio calculation of the photoionisation cross section of He has been performed 
for photon energies covering the non-resonant part of the spectrum from the ionisation 
threshold until about 300 eV. An analytical high-energy model tail has been introduced 
and the ah initio data were used in order to determine the single fit parameter. With 
this tail it became possible to predict the photoionisation spectrum for arbitrarily 
large photon energies within the underlying non-relativistic dipole approximation. 
Furthermore, the first-order correction to the dipole approximation could be estimated 
analytically with the aid of the model tail. 

Our theoretical results agree extremely well with the ones obtained by different 
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theoretical approaches in the low and the high energy parts by Venuti et al [B] and 
Ivanov and Kheifets [8], respectively. Therefore, we confirm the pronounced deviation 
between theoretical and experimental results noted in those earlier works. Particularly, 
at the photon-energy range around 50 eV that is very relevant to present-day FELs like 
FLASH or high-harmonic sources the relative deviation is unambiguously larger than 
the error estimates of the experiment of Samson et al jU |2] . Since we also demonstrated 
the validity of the dipole approximation at least for photon energies up to 200 eV where 
the previous comparisons were performed, this possible source of disagreement between 
theory and experiment is excluded. We thus conclude that the quality of theoretical 
calculations of the helium photoionisation cross section has reached a consistently higher 
level than the experiment. In view of the fundamental importance of the photoionisation 
cross section of helium we hope that the present work stimulates future experimental 
efforts to resolve the discrepancies between theory and experiment. Until this has been 
achieved, we propose to use the theoretical results, especially the ones of the present 
work that cover a large photon energy range in a consistent fashion, instead of the 
experimental ones as reference data for, e.g., calibrating new light sources, since they 
appear to be more accurate and reliable. 
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